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ABSTRACT 

Sampling  efficiencies  are  calculated  for  an  aspirated  parti- 
culate matter  sampling  probe  under  various  conditions  of  anlsoklnetlc 
flow.  A mathematical  model  developed  for  the  purpose  was  used  to  obtain 
results  for  a wide  range  of  particle  sizes  and  flow  velocities.  The 
results  can  be  used  to  predict  or  correct  sampling  errors  In  field  or 
laboratory  experiments.  Using  the  same  test  parameters  as  In  previous 
experimental  tests  by  other  workers,  sampling  efficiencies  were  cal- 
culated and  the  results  were  found  to  agree  favorably  with  the  results 
of  the  experiments. 
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NOTATION 


particle  concentration  In  the  sample,  g cm-3 
particle  concentration  In  the  free  stream,  g cm"3 
particle  diameter,  cm 

distance  from  the  Inlet  to  the  outlet  cross  section  of  the 
collection  tube,  cm 

thickness  of  the  collection  tube  wall  at  the  outlet  cross 
section,  cm 

length  of  coaxial  boundary  tube,  cm 
radial  co-ordinate  of  particle  position,  cm 
radius  of  coaxial  boundary  tube,  cm 
radius  of  collection  tube,  cm 

radial  co-ordinate  of  particle  position  far  upstream,  cm 

far  upstream  radius  of  the  stream  tube  that  impinges  on  the 
collection  tube  circumference,  cm 

time,  seconds 

radial  component  of  local  fluid  velocity,  cm  sec"1 
axial  component  of  local  fluid  velocity,  cm  sec"1 
fluid  velocity  in  collection  tube,  cm  sec"1 
fluid  velocity  at  boundary  tube  entrance,  cm  sec"1 
fluid  velocity  at  collection  tube  exit,  cm  sec"1 
fluid  velocity  at  boundary  tube  exit,  cm  sec"1 
free  stream  velocity,  cm  sec"1 

radial  component  of  local  particle  velocity,  cm  sec"1 

axial  component  of  local  particle  velocity,  cm  sec"1 

axial  co-ordinate  (origin  at  collection  tube  Inlet)  of  particle 
position,  cm 

axial  co-ordinate  of  particle  far  upstream,  cm 


NOTATION  (Cont'd) 


absolute  viscosity  of  fluid,  poise 
fluid  density,  g cm"3 
particle  density,  g cm"3 
stream  function,  cm3  sec"1 


The  following  are  dimensionless 


drag  coefficient  for  spheres 


dependent  variablessol ved  for  by  numerical  integration. 
They  represent  v^,  v^,  z and  r respectively 


collection  efficiency  of  sampling  tube 

thickness  of  collection  tube  wall,  h/rA 

grid  point  co-ordinates  In  the  radial  and  axial  directions 
respectively 

grid  point  co-ordinates  of  the  edge  of  the  collection  tube  Inlet 

axial  grid  point  co-ordinate  of  a particle  at  the  far  upstream 
position 

Inertia  parameter  of  particle 
radial  co-ordinate  of  particle,  r/rg 

radial  co-ordinate  of  particle  position  far  upstream,  rp  ^/r^ 

far  upstream  radius  of  the  stream  tube  that  impinges  on  the 
collection  tube  circumference,  r$  ^/r^ 

radial  co-ordinate  used  In  calculating  the  stream  function 
field,  r/rA 

spherical  particle  Reynolds  number  in  flow  In  the  proximity  of 
the  collection  tube 

spherical  particle  Reynolds  number  In  free  stream 


NOTATION  (Cont'd) 


up  radial  component  of  local  fluid  velocity,  du/dr 

uy  axial  component  of  local  fluid  velocity,  du/d z 

Vp  radial  component  of  local  particle  velocity,  dr/dT 

Vy  axial  component  of  local  particle  velocity,  dz/dt 

z axial  co-ordinate  (origin  at  collection  tube  Inlet)  of  particle, 

z/rB 

zQ  axial  co-ordinate  of  particle  far  upstream,  zQ/rB 

Z axial  co-ordinate  used  In  calculation  of  the  stream  function 

field,  z/rA 

a ratio  of  collection  tube  radius  to  boundary  tube  radius,  rg/rA 

8 length  of  coaxial  boundary  tube,  L/rA 

y distance  from  the  inlet  to  the  outlet  cross  section  of  the 

collection  tube,  0/rA 

6 distance  from  inlet  of  boundary  tube  to  inlet  of  collection 

tube,  B - y 

t time,  tUA/rg 

$ dimensionless  group  independent  of  particle  position,  ReQ2/K 

y stream  function,  4>/SUArA2 
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1.  INTRODUCTION 

The  collection  of  a representative  sample  of  finely  divided 
particulate  matter  from  still  or  moving  airstreams  is  required  where  the 
size  distribution,  mass  flow  rate,  concentration,  or  some  other  charac- 
teristic of  the  particulate-air  system  has  to  be  determined.  A sample 
of  particulate  matter  will  be  representative  only  if  the  particle  size 
distribution  and  content  in  the  sample  are  the  same  as  those  in  the 
ambient  air  at  the  point  of  sampling. 

Particulate  matter  sampling  instruments  used  at  DRES  consist 
essentially  of  a probe  or  tube  through  which  the  sample  is  drawn  and 
then  separated  from  the  air  stream  for  analysis. 

Such  a sampling  system  may  be  subject  to  three  distinct  types 
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of  error  (Vltols,  1964)  due  to: 

(1)  particles  failing  to  enter  the  sampling  probe  in 
representative  concentrations; 

(2)  particles  being  deposited  between  the  probe  mouth 
and  the  separation  location;  and 

(3)  particles  being  shattered,  aggregated,  or  incompletely 
retained  by  collection  devices. 

Particles  enter  the  sampling  probe  in  representative  concen- 
trations when  the  entrance  velocity  is  exactly  equal  to  the  velocity  of 
the  gas  being  sampled,  in  which  case  the  sampling  is  said  to  be  isokinetic. 
If  the  velocities  are  unequal  errors  of  type  (1)  occur,  in  which  case 
the  sampling  is  referred  to  as  anisokinetic. 

The  purpose  of  this  report  is  to  describe  a mathematical  model 
for  calculating  the  error  due  to  anisokineticity  for  a sampling  probe 
developed  and  used  at  ORES  (Fig.  1).  The  nominal  inside  diameter  of  this 
probe  is  3/4  inch.  The  inlet  end  is  sharp  edged  and  the  outside  surface 
is  sloped  so  that  the  tube  wall  thickness  increases  away  from  the  inlet 
with  a 1 in  12  slope.  The  length  of  the  tapered  section  is  2 3/4  inches. 
Results  are  provided  for  various  flow  conditions  so  that  the  user  of  the 
sampling  probe  can  determine  the  magnitude  of  sampling  errors  and  make 
corrections  when  anisokinetic  sampling  occurs. 

2.  DEFINITION  OF  THE  PROBLEM 

Due  to  inertial  and  drag  forces,  a particle  flowing  at  the  free 
stream  velocity  far  upstream  of  the  sampling  probe  will  not  necessarily 
follow  a stream  line  in  the  vicinity  of  the  probe,  where  radial  and  axial 
velocity  components  of  the  fluid  may  be  changing  markedly.  Thus,  if  the 
estimated  concentration  of  particles  in  the  free  stream  is  taken  as  the 
number  of  particles  collected  by  the  sampler  divided  by  the  volume  of  air 
passing  through  the  probe,  the  calculated  values  may  differ  markedly  from 
the  true  free  stream  values.  The  total  volume  of  air  passing  through  the 
probe  will  be  that  enclosed  by  a stream  tube  that  impinges  on  the  outer 
circumference  of  the  probe  opening.  When  the  velocity  at  the  tube  inlet. 
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Ug,  is  less  than  the  free  stream  velocity,  Uft  (Ug/UA  < 1),  particles  from 
outside  the  limiting  stream  tube  will  enter  the  sampler,  while  for  Ug/UA 
> 1 particles  originally  inside  the  limiting  stream  tube  will  pass  outside 
the  probe.  All  particles  of  the  same  diameter  that  are  collected  by  the 
probe  are  those  within  the  circular  envelope  generated  by  particles 
that  just  impinge  on  the  outer  circumference  of  the  probe.  Let  r be 

P,oo 

the  upstream  radius  of  the  limiting  particle  trajectory  envelope,  and 

be  the  upstream  radius  of  the  stream  tube  that  impinges  on  the  probe 
circumference.  Then  the  sampling  efficiency  is  given  by: 


where  CQ  is  the  particle  concentration  in  the  free  stream  and  C is  the 
particle  concentration  in  the  sample. 


The  collection  efficiency  of 


the  probe  is  given  by: 


(Eq.  2) 


where  rg  is  the  radius  of  the  probe  inlet.  The  problem  then  is  to  cal- 
culate r and  r$  ^ so  that  the  sampling  efficiency  can  be  found  from 
Eq.  1 and,  incidentally,  so  that  the  collection  efficiency  can  be  found 
from  Eq.  2. 


3.  EQUATIONS  OF  MOTION 


The  motion  of  an  individual  particle  has  been  shown  (Vitols,  1964 
and  Batchelor,  1956)  to  be  determined  by  the  following  ordinary  differential 
equations: 


dvF  CDRe^uF  ‘ 

dx  24  K 


(Eq.  3) 


dvj  C0Re(u^  - vj) 
dT  24  K 


(Eq.  4) 
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where 

Re  = ReQ[ (u^  - 

V)2  + % “ 

(Eq.  5) 

od2UA 

K = 

18urB 

particle  inertia  parameter 

(Eq.  6) 

U dp 

Re  = -2— 

0 V 

free  stream  Reynolds  number 

(Eq.  7) 

The  symbols  are  defined  in  the  notation  section  near  the  front 
of  this  report  and  the  basic  geometry  of  the  flow  system  is  illustrated 
in  Fig.  2. 


Several  assumptions  are  inherent  in  the  use  of  Eqs.  3 and  4 for 
calculating  the  collection  and  sampling  efficiencies  due  to  a stream  of 
particles,  including: 

(a)  uniform  particle  distribution; 

(b)  no  gravitational  or  electrostatic  forces  of  consequence; 

(c)  monodisperse  spherical  particles  with  diameter  very  small 
in  relation  to  the  inlet  diameter  of  the  probe;  and 

(d)  free  stream  flow  that  is  steady,  incompressible  and  irro- 
tational . 

The  drag  coefficient  is  a function  of  Reynolds  number  and  is 
available  in  the  form  of  definitive  empirical  equations  (Davies,  1945). 
These  equations  are  stated  as  follows: 


Re 

for  Re 
log10Re 

for  3 < 


= — 2.3363  x 10'4  (CnRe2)2  + 2.0154  x 10-6 

24  u 

6.9105  x 10-9  (CDRe2)4 

< 4 or  CpRe2  < 140 

« - 1.29536  + 9.86  x 10"1  (logioCpRe2)  - 4.6677 
(log10CDRe2)2  + 1.1235  x 10"3  (log10CDRe2)3 

Re  < 104  or  CpRe2  < 4.5  x 107 


(CDRe2)3  - 
(Eq.  8) 

x 10"2 

(Eq.  9) 
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AIR  FLOW  FIELD  EQUATIONS 


The  equations  of  fluid  velocity  were  derived  from  the  stream 
function  for  ideal  flow  over  and  through  the  collection  probe.  To  solve 
the  problem,  an  outer  boundary  was  used  around  the  collection  tube  in  the 
form  of  a coaxial  tube  of  radius  r^  (Fig.  3),  which  was  chosen  large 
enough  so  that  the  effect  of  the  boundary  tube  on  flow  in  the  proximity 
of  the  collection  tube  is  negligible.  The  collection  tube  was  inserted 
a distance  D into  the  downstream  end  of  the  boundary  tube.  Since  the 
flow  is  axisymmetric  only  a radial  plane  containing  both  tubes  has  to 
be  considered. 

The  fluid  enters  the  boundary  tube  with  steady  velocity  Uft, 
and  separates  into  a central  stream  with  velocity  Ug  at  the  entrance 
of  the  collection  tube,  and  an  annular  stream,  with  velocity  U^,  at 
the  downstream  end  of  the  boundary  tube.  The  axial  velocities  U^,  Ug 
and  Uc  are  uniform.  Also,  there  is  no  radial  flow  at  the  end  cross 
sections. 

The  boundary  conditions  on  the  flow  can  now  be  completely  speci- 
fied so  that  the  flow  field  can  be  obtained  by  solution  of  the  equation 
of  the  stream  function. 


The  axially  symmetric  stream  function  »p( r ,z ) (Batchelor,  1967) 


satisfies: 


iii.lM  + ifi,  o 

3r2  r 3r  3z2 


(Eq.  10) 


The  two  velocity  components  (Fig.  2)  are  given  by: 

,,  -1  ii 


1 li 


(Eq.  11) 


(Eq.  12) 


When  UA  and  Ug  are  specified  continuity  gives  Uc  as  follows: 
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>B  + 

-\  ^ 


(Eq.  13) 


where  a ■ 


(Eq.  14) 


and  h Is  the  thickness  of  the  collection  tube  wall  at  the  outlet  end 
(Fig.  3). 


For  uniform  velocity  profiles,  the  stream  function  Is  of  the 


* • Sur2  (Eq.  15) 

To  allow  for  greater  generality,  the  stream  function  and  the 
geometric  variables  were  restated  In  the  following  dimensionless  form: 


(Eq.  16) 


(Eq.  17) 


Z 


(Eq.  18) 


e - — 


(Eq.  19) 


(Eq.  20) 


(Eq.  21) 


The  boundary  values  for  the  stream  function  and  the  geometric 
configuration  In  terms  of  the  dimensionless  variables  are  shown  in  Fig.  4. 

The  axially  symmetric  stream  function  equation  (Eq.  10)  becomes: 


ifl  1 at 


3R2  R 3R  3R2 


(Eq.  22) 
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5.  DISCRETIZATION  SCHEME  FOR  THE  AIR  FLOW  FIELD 

The  equation  for  the  axially  symmetric  stream  function  (Eq.  22) 
Is  discretized  as  follows: 


aR2 


VlJ  , L.'i.-l.J  + vi,J-l 


21aR2 


AZ2 


(Eq.  23) 


where  1 and  j are  the  grid  point  numbers  in  the  R and  Z directions  res- 
pectively (Fig.  5).  Eq.  23  can  be  rearranged  to  give  a simple  equation 
by  choosing  a square  grid  so  that  AZ  and  AR  are  equal.  The  resulting 
equation,  which  Is  suitable  for  Gauss-Seldel  iteration  (Carnahan  et  al., 
1969),  Is  given  as  follows: 


, . . * . Tl+ij. vl-u 

1,j  4 81 


(Eq.  24) 


Eq.  24  Is  suitable  for  all  interior  points,  which  are  defined  as  points 
for  which  the  nearest  boundary  is  at  least  one  grid  size  unit  away.  A 
grid  can  be  chosen  such  that  all  points  not  right  on  the  boundary  are 
Interior  points,  with  the  exception  of  the  outside  boundary  of  the  col- 
lection tube  which  has  a slope  of  1/12.  To  handle  these  points,  which 
herein  are  called  special  boundary  points  (Fig.  6),  a Taylor  series  ex- 
pansion for  a point  near  the  boundary  was  used,  as  Is  generally  applied 
to  curved  boundaries  (Carnahan  et  al.,  1969).  The  resulting  finite 
difference  equation  is  suitable  for  these  points  Instead  of  Eq.  24. 


1.J 


2(a+l ) 


T1.J-1  + *1.j+l  + 


a(a+l)  a+1  1(l+a) 


(Eq.  25) 


where  a Is  the  distance,  in  terms  of  the  grid  size  unit,  from  the  grid 
point  to  the  sloping  boundary  at  point  A and  0 < a < 1. 


The  grid  size  was  chosen  so  that  the  thickness  of  the  collection 
tube  wall  at  the  downstream  end  corresponded  to  exactly  two  units.  Then 
the  length  of  the  tube  must  be  24  grid  units  since  the  slope  of  the 
outer  collection  tube  boundary  Is  1/12.  The  boundary  tube  radius 
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and  the  distance,  6,  (Fig.  4)  to  the  upstream  end  of  the  boundary  tube 
were  chosen  so  that  they  are  at  least  five  times  the  radius  of  the  Inlet 
to  the  collection  tube,  so  that  particles  travelling  in  the  flow  toward 
the  tube  inlet  will  behave  as  if  they  are  coming  from  far  upstream  and 
into  the  collection  tube  in  free  space.  Now,  to  satisfy  these  conditions, 
and  to  ensure  adequate  core  storage  space  in  the  DRES  IBM  1130  computer, 
while  maintaining  a fineness  of  grid  for  sufficient  accuracy,  the  col- 
lection tube  radius  was  chosen  to  be  12  units,  the  boundary  tube  radius 
60  units,  and  the  boundary  tube  length  88  units  (Fig.  5). 

The  stream  function  at  each  point  was  then  obtained  by  means 
of  Gauss-Seidel  iteration  (Carnahan  et  al.,  1969)  using  Eq.  24  and  Eq.  25. 
The  boundary  conditions  were  set  initially  and  held  fixed  throughout  the 
course  of  the  solution  (Fig.  4)  and,as  a starting  point  for  the  iteration, 
the  value  of  the  stream  function  at  all  points  not  right  on  the  boundary 
was  set  to  zero.  The  calculations  were  done  with  the  DRES  IBM  1130 
computer  by  means  of  a Fortran  program,  the  listing  of  which  is  shown 
in  Appendix  A. 

6.  SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

As  previously  stated  in  Section  2 of  this  report,  the  problem 

is  to  find  the  upstream  radii  r and  r so  that  the  sampling  and  col- 

P i00  s *°° 

lection  efficiencies  can  be  calculated.  In  the  same  dimensionless  form 

used  in  Eqs.  3 and  4 the  value  of  rp  m (notation)  was  found  by  an  iterative 

procedure  called  the  half  interval  method  (Carnahan  et  al.,  1969).  The 

value  of  r for  a critical  particle  was  estimated  far  upstream,  the 
P *°° 

path  followed  to  the  plane  of  the  collection  tube  opening,  and  the  miss 
distance  calculated.  Next,  the  half  interval  method,  previously  mentioned, 
was  applied  to  determine  a better  initial  estimate.  Then  the  path  was 
followed  again  to  the  plane  of  the  collection  tube  opening  for  another 
calculation  of  the  miss  distance.  This  process  was  repeated  several  times 
until  sufficient  accuracy  was  achieved.  The  initial  upstream  position  in 
a plane  perpendicular  to  the  flow  direction  was  located  far  enough  from 
the  collection  tube  opening  so  that  free  stream  conditions  prevailed.  A 
distance  of  five  target  radii  upstream  of  the  target  centre  was  considered 
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adequate  (Batchelor,  1956). 


The  path  of  an  Individual  particle  was  determined  step-by-step 
by  applying  a fourth  order  Runge-Kutta  method  (Carnahan  et  al.,  1969)  to 
the  equations  of  motion  (Eqs.  3 and  4).  The  values  of  Re  and  K In  these 
equations  were  easily  found  for  each  new  step  by  direct  substitution  of 
previously  determined  values  into  Eqs.  5,  6 and  7,  but  the  value  of 
CpRe  In  Eqs.  3 and  4 had  to  be  calculated  In  each  step  by  numerical  solu- 
tion of  the  definitive  empirical  equations  (Eqs.  8 and  9).  This  was  done 
using  Newton's  method  (Carnahan  et  al.,  1969)  for  finding  the  zero  of  a 
function.  The  values  of  u^  and  uj  were  calculated  in  each  step  from  the 
stream  function  field  as  follows: 

uL  - T1J-1  ' 
r 4 ( i - 1 ) (AR)2 


(Eq.  26) 


jjl,  !L*1xL1 
z 


4(1-1 )(AR)2 


(Eq.  27) 


where  i and  j define  the  grid  point  of  the  particle  position.  Since  the 
radius  of  the  collection  tube  was  chosen  to  be  12  grid  units  these  are 
given  by: 


1 = 1 + 12r 
j - jQ  + 12(1  - z0) 


(Eq.  28) 
(Eq.  29) 


where  j and  iQ  are  the  starting  point  values  of  j and  z.  The  values  of 

i and  j obtained  from  Eqs.  28  and  29  were  rounded  off  to  the  nearest 

lower  Integer  value  in  each  calculation.  The  value  of  r,  was  obtained 

s.» 

directly  from  the  stream  function  by: 


calculated  at  the  lowest  value  of  1 satisfying 
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where  ig  and  jg  define  the  grid  point  at  the  edge  of  the  collection  tube 

inlet.  The  calculations  to  obtain  the  solutions  were  done  with  the  DRES 

IBM  1130  computer  by  means  of  a Fortran  program,  the  listing  of  which  is 

shown  in  Appendix  B.  The  sampling  and  collection  efficiencies  given  by 

Eqs.  1 ana  2 were  also  obtained  by  this  program  after  the  values  of 

? and  F had  been  calculated. 
p,°°  s ,°° 

7.  RESULTS 

A sample  stream  function  field  for  one  set  of  input  data  is 
shown  after  its  associated  computer  program  (Appendix  A),  and  a sample 
calculation  of  the  sampling  and  collection  efficiencies  is  shown  after 
their  associated  computer  program  (Appendix  B).  Using  these  two  computer 
programs  many  more  calculations  were  made  to  produce  the  graphical  results 
shown  in  Figs.  7 to  14.  These  results  are  described  in  greater  detail  as 
follows. 

Many  experimenters  have  measured  sampling  efficiencies  over  the 
last  sixty-five  years  (Vitols,  1964).  An  experimental  study  using  zinc 
sphere  test  dust  in  a wind  tunnel  was  done  for  sampling  tubes  of  0.65  to 
1.90  cm  diameter  (Badzioch,  1959).  The  tubes  were  blunt-edged  with  a wall 
thickness  of  0.6  mm.  The  results  for  some  of  these  experiments  are  shown 
in  Fig.  7.  Using  the  same  test  parameters  as  in  these  experiments,  samp- 
ling efficiencies  were  calculated  for  various  velocity  ratios,  Ug/U^,  and 
the  results  plotted  in  Fig.  7 assuming  that  Ug/U^,  of  the  mathematical 
model,  is  equivalent  to  U/Uq  in  the  experiments  of  Badzioch.  The  results 
agree  favorably  even  though  the  tube  shapes  were  slightly  different.  How- 
ever, calculations  for  various  input  parameters,  keeping  Ug/U^  = 1,  gave 
sampling  efficiencies  between  0.9900  and  1.0000,  indicating  that  the 
tapered  outside  wall  of  the  collection  tube  does  not  disturb  the  flow  so 
as  to  cause  deviations  for  isokinetic  velocity  conditions.  Therefore, 
the  two  tubes  can,  for  practical  purposes,  be  considered  equivalent. 

The  sampling  and  collection  efficiencies  are  functions  of  two 
dimensionless  groups,  the  inertial  parameter,  K,  and  the  free  stream 
particle  Reynolds  number.  Re  . A new  dimensionless  group: 
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(Eq.  32) 


independent  of  particle  size  can  be  introduced  (Friedlander,  1977). 
According  to  the  rules  of  dimensional  analysis  this  is  permissible,  but 
the  efficiencies  are  still  determined  by  two  groups  chosen  to  be  K and  $. 
The  calculated  sampling  efficiency  is  plotted  against  the  inertia  para- 
meter in  Figs.  8 and  9 for  sampling  rates  of  5 and  40  litres  per  minute 
respectively.  Curves  are  plotted  for  values  of  Ug/U^  of  1/3,  1,  3 and 
10  in  each  of  the  two  figures.  Note  that  equal  values  of  Ug/Uft  correspond 
to  a different  value  of  $ in  Figs.  8 and  9.  This  is  due  to  two  different 
sampling  rates  because  sampling  rate  and  the  corresponding  free  stream 
velocity  were  the  only  dimensional  parameters  which  were  changed.  The 
calculated  collection  efficiencies  were  plotted  in  Figs.  10  and  11  for 
the  same  velocity  ratios  and  sampling  rates  as  the  sampling  efficiencies 
shown  in  Figs.  8 and  9. 

To  show  the  effect  of  varying  the  flow  velocity  in  the  sampling 
tube  while  maintaining  a constant  velocity  ratio,  the  sampling  efficiency 
was  plotted  against  the  inertia  parameter  in  Fig.  12.  Curves  were 
plotted  for  the  two  sampling  rates  of  5 and  40  litres  per  minute  as  before. 
Similar  curves  are  shown  for  the  collection  efficiencies  in  Fig.  13. 

To  illustrate  directly  the  drastic  effect  of  varying  free  stream 
velocity  on  the  sampling  efficiency  for  constant  sampling  velocity  and 
particle  size  the  sampling  efficiency  was  plotted  against  the  velocity 
ratio,  Ug/U^,  for  a particle  size  of  100  pm  in  Fig.  14. 

8.  DISCUSSION 

As  can  be  seen  in  Fig.  7,  where  the  experiments  of  Badzioch 
are  compared  to  the  results  obtained  by  means  of  the  present  mathematical 
model,  the  effect  of  anisokineticity  produces  sampling  errors  of  greater 
than  100  percent  even  when  the  sampling  velocity  differs  from  the  free 
stream  velocity  by  less  than  a factor  of  three.  These  results  are  for 
23  micron  spherical  particles.  Smaller  particles  give  rise  to  smaller 
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sampling  error.  The  exact  errors  for  very  small  particles  cannot  be 
determined  by  the  present  model  because  the  computing  errors  Increase 
with  decreasing  particle  size  for  very  small  particles.  Experience  shows 
that.  In  addition  to  the  requirement  of  a small  grid  size  for  calculating 
the  stream  function  field,  reliable  results  can  only  be  achieved  If  the 
values  of  the  time  step  Increment  dt  Is  at  least  less  than  the  Inertia 
parameter  K.  So,  as  particle  size  decreases,  and  correspondingly  K de- 
creases, a larger  number  of  calculations  are  required  which.  In  Itself, 
can  give  rise  to  computing  errors.  However,  the  results  did  Indicate 
that  the  sampling  errors  for  particles  In  the  order  of  10  microns  can  be 
at  least  50  percent.  For  example,  for  a sampling  rate  of  5 litres  per 
minute  and  a sampling  velocity  which  Is  a third  of  the  free  stream  velo- 
city, Fig.  8 shows  a value  of  the  sample  concentration  which  Is  1.5  times 
the  free  stream  concentration  when  the  Inertia  parameter  Is  0.029.  This 
corresponds  to  a particle  size  of  10  microns  when  the  particle  density  is 
1 g cm-3  and  5 microns  when  the  particle  density  is  4 g cm"3.  The  density 
of  zinc  sulfide,  which  Is  sometimes  used  as  phosphorescent  trace  material 
with  particle  sizes  of  about  5 microns,  is  4.1  g cm"3  (Sehmel,  1973). 

Although  the  model  cannot  be  used  for  calculations  using  completely 
still  air,  very  low  free  stream  velocities  can  be  treated.  For  a sampling 
rate  of  5 litres  per  minute  and  a velocity  ratio  Ug/U^  * 10,  as  used  In 
one  of  the  curves  in  Fig.  8,  the  corresponding  free  stream  velocity  Is 
3 cm  sec"1.  This  is  about  the  wind  velocity  expected  In  an  enclosure  such 
as  an  ordinary,  normally  ventilated,  room. 

Another  point  to  note,  when  considering  sampling  errors.  Is  that 
even  when  the  velocity  ratio  and  particle  size  are  kept  constant,  the 
sampling  error  varies  somewhat  with  the  absolute  value  of  the  sampling 
tube  velocity.  This  Is  indicated  by  the  difference  in  the  two  curves 
(Fig.  12)  for  sampling  rates  of  5 and  40  litres  per  minute. 

The  drastic  effect  on  sampling  efficiency  of  varying  the  free 
stream  velocity  while  maintaining  a constant  sampling  velocity  Is  illus- 
trated In  Fig.  14  where  the  calculated  results  for  a particle  size  of  100 
microns  are  plotted  directly. 
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9.  CONCLUSIONS 

The  effect  of  anlsoklnetlclty  on  sampling  with  the  DRES  sampling 
probe  used  In  field  conditions  Is  large  enough  so  that  drastic  errors  In 
sampling  can  occur.  The  results  obtained  by  the  mathematical  model  des- 
cribed herein  can  be  used  to  estimate  the  size  of  these  sampling  errors. 
They  can  also  be  used  to  correct  measured  samples  If  the  wind  velocity 
and  sampling  rate  are  also  measured,  and  the  particle  size  and  density 
are  known.  Long  sampling  periods  are  best,  so  that  effects  due  to  wind 
speed  changes  which  occur  over  short  time  periods  are  averaged  out. 
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Figure  1:  Collection  Tube  for  Sampling  Particulates  in  Air 
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FIGURE  2:  CO-ORDINATE  SYSTEM  FOR  AXIAL  FLOW 

IN  THE  PROXIMITY  OF  THE  COLLECTION  TUBE 
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FIGURE  12=  EFFECT  OF  FLOW  VELOCITY  IN  SAMPLING  TUBE  ON  SAMPLING  EFFICIENCY 

FOR  FIXED  VELOCITY  RATIO  (±J»-  = ^ ) 
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FIGURE  14:  EFFECT  OF  FREE  STREAM  VEUOCITY  ON  SAMPLING  EFFICIENCY 

FOR  ONE  PARTICLE  SIZE  AND  CONSTANT  LOW  VELOCITY  IN  SAMPLING  TUBE 
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"D"  In  Figures  7 to  14  denotes  the  diameter  of  the  collection 
tube.  In  terms  of  the  symbols  defined  under  Notation  this 
should  read  2rg. 

Figure  14  caption,  second  line,  should  read  "constant  flow 
velocity". 


